


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1963 


An analysis of the properties of 
two-dimensional incompressible fluid flow in 
the mixing chamber of a constant area ejector 


Kiefer, Martin D. 


Monterey, California: U.S. Naval Postgraduate School 
http://ndl.handle.net/10945/11630 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
(8 DUDLEY research materials and institutional publications created by the NPS community. 
«ist sae Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS'‘s first 


INN KNOX appointed — and published -- scholarly author. 

| LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


NPS ARCHIVE 


1963 
KIEFER, M. 





AN ANALYSIS OF THE PROPERTIES 
OF TWO-DIMENSIONAL INCOMPRESSIBLE 
FLUID FLOW IN THE MIXING CHAMBER 

OF A CONSTANT AREA EJECTOR 


MARTIN D. KIEFER 


LIBRARY 
US. NAVAL POSTGRADUATE SCHOOL 
MONTEREY CALIFORNIA 


DLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOO! 
MONTEREY CA 93943-5191 

















AN ANALYSIS OF THE PROPBRRTIS 
OF TWO-DIMENSTONAL INCOMPRESSIBLE FLUID FLOW 


IN THE MIXING CHAMBER OF A CONSTANT? ARHKA EJECTOR 


Martin D, hiefer 





AN ANALYSIS OF THE PROPERTIES 
OF TWO-DIMENSIONAL INCOMPRESSIBLE FLUID FLOW 
IN THE MIXING CHAMBER OF A CONSTANT AREA EJECTOR 


by 
Martin D, i one 


Lieutenant, United States Navy 


submitted in partial fulfillment of 
the requirements for the degree of 


MASTER OF SCIENCE 
ENG 
MECHANICAL ENGINEERING 


United States Naval Postgraduate School 
Monterey, California 


19 6 3 








| LIBRARY | 
U.S. NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA 


AN ANALYSIS OF THE PROPERTIES 
OF TWO-DIMENSIONAL INCOMPRESSIBLE FLUID FLOW 
IN THE MIXING CHAMBER OF A CONSTANT AREA EJECTOR 


by 
Martin D. Kiefer 


This work is accepted as fulfilling 
the thesis requirements for the degree of 
MASTER OF SCIENCE 
IN 
MECHANICAL ENGINEERING 
from the 


United States Naval Postgraduate School 





ABSTRACT 

The properties of an incompressible fluid in the mixing 
chamber of a constant area cylindrical ejcctor are analyzed 
in a twoedimensional form. The flow field for velocity and 
temperature is assumed to be made up of regions in which 
uniform flow of the orimary and secondary fluids exist, and 
which diminish with axial distance from the ejector entrance, 
The velocity and temperature distributions in the mixed 
region are assumed to have cosine shaped profiles, the con- 
patible solutions are given an axial distribution by assuming 
& parabolic spread rate for the secondary fluid fet boundary. 

Results are venerated in the form of a nonedimensional 
velocity and pressure, 2 normalized temperature and a omen- 
twn Factor, as functions of axial distance, Area ratios 
studied include 10031, 9:1 and 2,263:1. Velocity ratios 
studied include 50:1, 10:1, 5:1 and 1.5:1. The Fortran vro- 
grams employed to generate compatible solutions, titled 


RJECTMIXN I and SJECTMIX II, are included, 
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1. Introduction 

The majority of analyses of the physical properties of 
fluids in ejectors have used a one-dimensional approach in 
which uniform properties are assumed at the entrance to the 
mixing chamber which is of sufficient length for the pro- 
perties to again become uniform at the exit.t For example, 
consider a primary fluid flowing along the central axis of a 
cylindrical ejector and a secondary fluid flowing parallel to 
and completely surrounding it, as shown in Fig. 1. These 
fluids are separated by a circular boundary. At the entrance 
to the mixing chamber the boundary between the two fluids is 
removed and, at this initial cross-section, both the primary 
and the secondary fluids have uniform but different velo- 
cities, As the primary and secondary fluids flow axially 
through the mixing chamber, the fluids mix and the physical 
properties of each one of them are affected by the physical 
properties of the adjacent fluid. 

It is the purpose of this investigation, by use of a 
high speed digital computer and Fortran Programming Language, 
to predict what the properties in the partially mixed, non- 

ls Pai, Fluid Dynamics of Jets, D. Van Nostrand Co., 
Inc., New York, 195i 

eR, A, Smith, "Theory and Design of Simple Ejectors", 


=e Aspects of Fluid Flow, Edward Arnold and Co., London, 
Lio 
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uniform sections of the stream will be. Consideration of 
the fluid flow mechanisms, such as the shear stress within 
the fluids, leads to the concept of a profile shape for the 
properties at a given cross-section. A typical profile shape 
for velocity and temperature is shown in Fig. 2. 

The first problem considered is one in which the velo- 
cities of the primary and secondary fluids are different but 
their temperatures and densities are the same. The Fortran 
Program developed to solve this problem is called EJECTMIX I. 

The second problem, considered in the program BJECT- 
MIX II, studies the variations in the velocities and the 


temperatures of the two fluids, 


3p, D,. McCracken, A Guide to FORTRAN Programming, 
John Wiley and Sons Inc., New York, 1961 





2. Ejector Analysis 

In order to analyse the behavior of the physical proper- 
ties in the mixing chamber of an ejector there are several 
possible analytic approaches. 

The first possible approach is a one-dimensional over- 
all analysis in which uniform properties are assumed at the 
entrance and uniform properties of the completely mixed 
fluids are assumed at the exit of the mixing chamber. With 
this analysis, the region in which mixing of the primary and 
secondary fluids takes place is not considered. The only re- 
gions where the physical properties are studied are at the 
entrance to the mixing chamber and after the point where uni- 
form flow exists. 

A second possible approach is that of an internal analy- 
sis which is based on the actual fluid mechanisms. This ap- 
proach requires the examination of an element of fluid and 
the forces which are acting upon it. The difficulty here 
lies in the inability to acquire a usable analytic expression 
for the shear stress acting upon the element of fluid, 

A third approach, the one used in this thesis, is the 
generation of compatible solutions for the physical proper- 
ties of the fluids at representative axial cross-sections by 
means of consideration of Conservation of Mass, Momentum, 
Bernoulli's Equation and the Energy Equation. In order to 
apply the above equations to the fluid in the mixing chamber, 


certain basic assumptions as to the nature of the flow in the 
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mixing chamber must be made, 

First of all, the flow in the mixing chamber is assumed 
to be divided into regions which are shown in Fig. 3a and 
Fig. 3c. These regions are defined by the central core 
boundary and the expanding jet boundary. There are similar 
velocity and temperature profiles which consist of uniform 
portions at the central axis and the periphery of the mixing 
chamber, which are connected by a cosine shaped profile. The 
assumed profile shapes are shown in Fig. 3b and Fig. 3d. 
Assuming profile shapes, negligible boundary effects, con- 
stant pressure across any given cross-section, and constant 
temperature until a mixing region is reached, the previously 
mentioned equations can then be used, 

The first equation considered is the Equation of Contin- 
uity 

me f's V aA =Constant 
O 
For a typical profile shape, as shown in Fig. 3e, the 
Equation of Continuity becomes 
r 1 r 
m= 2p [Vp v dr +f VYave* av coo (Fe | i: dr +f “Vs yr dr 
Oo ag MH -Ya YM, 
Integrating and evaluating the integral gives 
m = Tp Vp r_-+Vave Pye-r a4 - hav [peal +Vs Pen] 


Dividing through the entire equation by 1pVsra’, the 
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equation takes on the non-dimensional form 





The next consideration is that of Bernoulli's "Equation 
alone a streamline in the uniform velocity portions of the 


profile. Bernoulli's Fquation states that 


pve = Constant 
PrP 2 


In non-dimensional form the equation hecomes 


p Pole -Constant 
2 pV s' Ta? 


The constant associated with Bernoulli's Equation in the 


primary core is found by evaluating Bernoulli's Equation in 
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the non-dimensional form along the axis at the antrance to 
the mixing chamber, Similarly, the constant associated with 
the secondary cone is determined at the entrance to the mix- 
ing chamber, These values will be used in conjunction with 
tone ilomentum “quation, 

Considerine the same profile which was used in the 


Equation of Continuity and solving for the Homentum rate 


ifomentum rate -l p Ve aA 
over the entire cross-section the Momentum Hate Squation be- 
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. a 2 
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Integrating and evaluating the integral 
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The Momentum Equation is then applied to the ejector in 
the axial direction. Neglecting shear stress at the ejector 
walls, 


Ap nr? = A(Momentum Rate) 


Defining a non-dimensional Momentum Rate by dividing by 


2 z : 
ees. ri» we obtain, 
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Ap, ts 2 AN 
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Definine a non-dimensional pressure difference, An*. as 
Rap) 3 3 


A pt. ae Ae Bray. 
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where Op =px-Pxe, and AM =M, —My eos 

The final equation used is the Hnergy Equation in which 
the neat transferred to the surroundings is assumed to be 
zero, The difference in geo-potential energy is negligible 
and, lastly, the chance in kinetic enerry is negligible in 
comparison with the Snthalpy Rate. Thus the “nergy Squation 
becomes 


it = fail = f nar ae Men [TV rdr =Constant 








Using the same profile shapes as before, the Enthalpy 


Rate for a typical axial cross-section becomes 
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was solved by numerical integration. The interval between 


van and The was divided into 100 increments and summed, 


The other integrals which make up the total value of the 


Enthalpy Rate are integrated and evaluated giving 
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In order to place this in non-dimensional form, the 
equation was divided by Tp ¢, Mig a Thus the integral 
o “oe 


for numerical integration becomes 


yt [V 
vex, ay Cos/ P7Pax dT, t alr 
t, Ae Vs ps On| Ea 1) |[Pavest os( ate aE r | 


Non-dimensionalizing the evaluated portion of the 


Enthalpy Rate gives 
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3. Method 

The equations employed in EJECTMIX I and BJBCTMIX II 
were used in the following manner. At a given axial position, 
a core boundary (r,) is assumed. Primary anc secondary velo- 
cities Ue and V,) are assumed which will be compatible with 
the Bernoulli Equation for the entrance to the mixing cham- 
ber, A jet boundary (xr,) is assumed and the conditions are 
checked for Conservation of Mass by means of the Continuity 
“quation. The primary and secondary velocities are corrected 
by an iterative process until a compatible solution is found. 

The velocities found by the Continuity tauation sre then 
used to solve for the pressure at the cross-section under 
consideration. Pressure changes determined from Bernoulli's 
aquation along streamlines in the uniform flow velocity re- 
sions are compared with the pressure change found from the 
overall Nomentum Rate. If the pressure is found to be incom- 
patible, the jet boundary (r,) is correctcd and once again 
the Equation of Continuity is used in order to determine 
compatible central and peripheral velocities. The values 
round by the Continuity Equation are then vsed to find a new 
value of the Momentum Rate, which, in turn, determines a new 
pressure change which can be used in Bernoulli's Equation to 
check the velocities found by Continuity. This process re- 
peats itself until the core and jet boundaries ie and r,), 
the central and peripheral velocities we and V.) and the 


pressure change are all compatible in the three bisic equa- 
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Fig.4 Block Diagram for EJECTMIX | 
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tions used, Then the process is repeated for the next cross- 
section which is displaced axially in the direction of fluid 
flow, 

_ After several cross-sections have been investigated, the 
core boundary (r,) or the jet boundary Oar) may reach their 
limiting value, that is, the central core boundary (r) may 
go to zero or the peripheral jet boundary (r,) may reach the 
mixing; chamber boundary (ro). 

Considering the case of either one of the boundaries 
reaching its limiting value while the other boundary remains 
within its limit, the scheme for finding the compatible solu- 
tion remains exactly the same as that used when both bound- 
aries were within limits, with the exception of the role of 
Bernoulli's “auation, Since the flow field will no longer 
be uniform where the boundary has reached its limit, the 
streamline along which Bernoulli's Iquation was assumed to 
have a constant value will no longer exist, With this in 
mind, Bernoulli's Equation will only be taken into account in 
the uniform flow portion which has not reached its Limiting 
value, 

When both boundaries have reached their respective limit- 
ing values, no additional restrictions are required other 
than the overall Momentum Rate. Compatible solutions are 
then found using fixed increments of the central velocity as 


the criteria for the particular cross-section. 


ay 


JECTHMIA II ts basically the same prorram as SJSCTMIA I 


15 











——————_—_—_—— 


EEE 


START 





iy 
(Total Energy), (Total Energy), 






(Total Energy), 


olutio 
Correc 
Change Rat 


Fig.5 Block Diagram for EJECTMIX II 


Yes Aolution 
“orrec 


Change Tp Change Ts 


Next Section 





16 





with additional restrictions placed upon the flow field, 
These additional restrictions take the form of variations 
in the temperatures of the primary and secondary fluids. The 
temperature boundaries are assumed to spread more rapidly 
than the velocity region boundaries, This has been verified 
by experiments on jets of gases flowing in ejectors,4+ Also, 
the temperature is assumed to be a constant until the temp- 
erature boundaries reach their limiting values. The temper- 
ature boundaries are assumed and the value of the energy equa- 
tion at the cross-section in question is compared with that 
of the initial cross-section. The central temperature core 
boundary (r,) is corrected until a compatible solution is 
obtained, 

In order to determine the axial position on which each 
of the representative compatible solutions lie, a spread rate 


E inp] is assumed for the peripheral jet boundary. When 
dx 


this boundary reaches its limiting value, it is assumed that 
the pressure will continue varying as a smooth curve until 


it reaches its limiting value in the mixing chamber, 


My, Kuchemann and J. Weber, Aerodynamics of Propulsion, 
McGraw-Hill Book Co., Inc., New York, 1953 
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i, Results 

The data which was obtained using EJECTMIX I and 
EJECTMIX II is presented in eraphical form in Appendix 2 and 
in tabular form in Appendix 3, The properties obtained are 
presented in groups according to the ratio of the radius of 
the peripheral boundary to the central jet radius, The 
values of this given ratio are 1.5, 3.0 and 10.0, 

These graphs may be entered with the known physical 
dimensions of the ejector mixing chamber and the velocity 
ratio, The ratio of the initial perinheral radius to the 
central jet radius of ten to one gives an actual area ratio 
of 100 to one and thus approaches the condition whereby the 
peripheral fluid actually would be unbounded, 

The graphical solutions shown are those in which the 
fluid was assumed to return to a uniform flow after it had 
passed down an axial distance equal to ten times the peri- 
pheral diameter. The data and graphs are presented as func- 
tions of axial distance divided by initial central jet radius 
(X/r). Tabulated data shown in Appendix 3 is the computer 
output which was used to produce these graphs. 

In the course of finding compatible solutions, a temper- 
ature difference of 100 degrees was used. 8y varying only 
the temperature, the solutions obtained for temperature were 
directly proportional to each other, With this in mind, the 
temperature difference for any case can be determined by 


using a proportionality constant and a 100 degree temperature 
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difference, Temperatures are measured above a reference 
temperature which is the secondary fluid temperature at the 
initial cross-section, 

The pressure term as used in Bernoulli's Equation is not 
significant in its absolute value since it cancels out of the 
equation when used in these programs, Its significance lies 
in the change of the non-dimensional pressure term, With 
each specific case an initial non-dimensional pressure can 
be determined, This non-dimensional pressure will be a func- 
tion of the absolute velocities, the pressure and the fluid 
flowing in the ejector, 

Since the initial peripheral velocity ratio is equal to 
one, a change in the non-dimensional pressure of one, before 
the velocity jet boundary reaches its limiting value, leads 
to a zero velocity ratio through the Bernoulli “quation. Any 
increase in the non-dimensional pressure beyond the value of 
one leads to compatible solutions in which the peripheral 
velocity becomes negative. This phenomenon is dependent 
upon the velocity ratio and the cross-sectional area of the 
central and peripheral jets. 

Negative peripheral velocity ratios first appeared in 
the range of the values studied at an initial central velo- 
city ratio of ten and an area ratio less than 100 to one, 

The flow pattern which exists when negative peripheral velo- 
city ratios are generated in the compatible solution is 


shown in Fig. 6, 
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Fig. 6 Profiles and Flow Pattern 
with Negative Peripheral Velocities 
Throughout most of the program accuracy of one-half of 
ons percent is required in order for the solution to be con- 
sidered acceptable, This criterion fails during the transi- 
tion of the peripheral velocity ratio term through its zero 
values. For this portion of the program accuracy of plus or 
minus one one-hundredth for the non-dimensional velocity term 
is required, This change in the criterion ensures that a 
finite and acceptable error can be established for any value 
of the peripheral velocity ratio. 
The final piece of data presented by the program is the 
Momentum Factor. This term gives that multiple of the final 


momentum that exists at any point. 
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5. Conclusions 

The programs presented here represent the first step in 
a series of programs which would eventually describe the 
physical properties of any fluid flowing in an ejector, This 
logical. continuation would require that work be done on the 
compressible fluid problem at both subsonic and supersonic 
velocities. 

The programs contained in this work solve the ejector 
problem for two incompressible fluids of the same density 
Whose physical properties are initially different. The data 


obtained is presented here in graphical and tabular form, 
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APPENDIX I 
Appendix I contains a listing, in FORTRAN language, of 
the programs EJECTMIX I and HJECTMIX II. 
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APPENDIX II 
Appendix II is a graphical presentation of the data 
obtained by programs EJECTMIX I and EJHCTMIX II, Various 
physical characteristics are plotted against the axial dis- 
tance, The axial distance is expressed as a multiple of the 
central jet radius, The various curves shown represent data 


obtained using different initial central velocity ratios, 
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Multiple of Final Momentum Rate vs, Axial Distance 
for an Ejector Area Ratio of 100:1 
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APPENDIX III 
Appendix III contains edited data from which the graphs 
of Appendix II were drawn, The headings are the sane as 


defined in the section of Symbols and Abbreviations, 
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